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O ■ Abstract 

, ' The identihcation of the X(3872) as a loosely-bound charm- meson molecule allows it to be 

j^ ■ described by an effective field theory, called XEFT, for the D*D, DD* and DDir sector of QCD 

at energies small compared to the pion mass. We point out that this effective field theory can be 
extended to the sector that includes an additional pion and used to calculate cross sections for the 
scattering of a pion and the X(3872). If the collision energy is much smaller than the pion mass, 
the cross sections are completely calculable at leading order in terms of the masses and widths of 
r^ • the charm mesons, pion masses, and the binding energy of the X(3872). We carry out an explicit 

calculation of the cross section for the breakup of the X(3872) into D*^D*^ by the scattering of a 



p. 



^H, very low energy vr^ 



> 

OO 

oo 

in 

o 
o 



X 



PACS numbers: 14.40.Rt, 13.75.Lb, ll.30.Rd 

Keywords: effective field theory, pion scattering, exotic mesons 



I. INTRODUCTION 

The X(3872), which was discovered by the Belle Collaboration in 2003 [1], is a truly 
remarkable hadron. There is increasingly compelling experimental and theoretical evidence 
that it is a loosely-bound charm-meson molecule whose particle content is 

X = ^ (D*0D° + D°D*°) (1) 

v2 

and whose constituents have a mean separation much larger than normal hadrons. Its 
J^*" quantum numbers are 1^^. Its binding energy, Ex, relative to the D*^D^ threshold 
is extremely small, less than 1 MeV. As a consequence, it has universal properties that 
are determined by Ex or, equivalently, by the large D*^D^ scattering length in the even 
charge-conjugation channel [2]. 

Fleming, Kusunoki, Mehen, and van Kolck have constructed an effective field theory 
called XEFT that can be used to calculate corrections to the universal properties of the 
X(3872) systematically [3]. In this effective field theory, the elementary constituents are 
the neutral charm mesons -D°, -D*", /)", and Z)*° and the neutral pion 7r°. There are two 
types of interactions: contact interactions between the pairs of charm mesons D*^D^ and 
D^D*^ and the pion transitions D*° <-)■ D°7r° and Z)*° ^ D^ti'^. XEFT was designed to 
describe systems consisting of D*^D^, D'^D*^, and D^D^tt^ with total energy very close to the 
D*0£)0 threshold, including the X(3872). In Ref. [3], the authors presented power-counting 
arguments that guarantee that the pion transitions can be treated perturbatively. They 
calculated the decay rate of the X(3872) into D^D^tt^ to next-to-leading order in XEFT. 
Fleming and Mehen have applied XEFT at leading order to decays of the X(3872) into the 
P-wave charmonium states Xcj and one or two pions [4]. They factored the amplitude for 
the decay into a long-distance XEFT matrix element and a short-distance factor that can 
be calculated using heavy-hadron chiral perturbation theory (HHxPT). 

The original formulation of XEFT has a rather limited domain of validity. One limitation 
is that it does not describe charged charm mesons. This limits its domain of applicability to 
the energy region within a few MeV of the D*^D^ threshold, because the D*^D~ threshold 
is higher only by about 8 MeV. This limitation is easily removed by generalizing XEFT to 
include the charged charm mesons D^ ^ D*^, D^ , and D*~ and the charged pions tt"^ and 
7i~ . This generalization extends the domain of applicability of XEFT to all energies relative 
to the D*^D^ threshold that are small compared to the D* — D mass difference, which is 
approximately equal to the pion mass 771^^ ^ 135 MeV. If the P-wave charmonium state 
Xci(2P), whose quantum numbers are also 1"'""'", lies in this region, it may be necessary to 
also include it as an explicit degree of freedom in order to carry out accurate calculations 
beyond leading order. 

Canham, Hammer, and Springer have obtained universal results for scattering of charm 
mesons with the X(3872) that depend only on the D*^D^ scattering length [5]. They 
calculated the S-wave phase shifts for D^X scattering and for D*^X scattering by solving 
the three-body problem for the charmed mesons. The D^X scattering length and the D*^X 
scattering length are both about an order of magnitude larger than the D*^D^ scattering 
length. These universal results can be regarded as the predictions of XEFT at zeroth order 
in the pion transitions. Thus XEFT can be applied to systems consisting of three charm 
mesons with energy close to the appropriate threshold. 

In this paper, we point out that XEFT can also be applied to systems consisting of D*D*, 
D*D7i, DD*7T, and DDnn with total energy close to the D*D* threshold. The scattering 



processes whose cross sections are calculable using XEFT include nX elastic scattering and 
ttX — )■ D*D* at collision energy much smaller than m,r- In Section II, we summarize the 
case for the X(3872) as a loosely-bound charm-meson molecule and we describe some of its 
universal properties. In Section III, we calculate the cross section for tt+X — )• D*^D*^. In 
Section IV, we discuss elastic tt+X scattering. Our results are summarized in Section V. 

II. THE X(3872) 

We begin by summarizing the case for the X(3872) as a loosely-bound charm- meson 
molecule whose particle content is given in Eq. (1). The only experimental information that 
is necessary to make this identification is the determination of its quantum numbers and the 
measurements of its mass. The quantum numbers of the X(3872) are l"*"^, which follows 
from 

• the observation of its decay into J/ip'~f, which implies that it is even under charge 
conjugation [6], 

• analyses of the momentum distributions from its decay into J/ipTr^n' , which imply 
that its spin and parity are 1"*" or 2" [7], 

• either the observation of its decays into D^D^tt^ [8], which disfavors 2~ because of 
angular momentum suppression, or the observation of its decay into ip{2S) 7 [9], which 
disfavors 2~ because of multipole suppression. 

Measurements of the mass of the X(3872) in the J/ipTf^TC^ decay channel [10] imply that 
its energy relative to the D*^D^ threshold is 

Mx - (Mzj.o + Mzjo) = -0.42 ± 0.39 MeV. (2) 

The mass of the X(3872) has also been measured in the D^D^n^ decay channel [8] and in the 
£)*0£)0 dQQoy channel [11]. Measurements in the D^D^tt^ channel are biased towards larger 
values by a contribution from a threshold enhancement of D*^D^ above the threshold [12]. 
Measurements in the D*^D^ channel are further biased towards larger values by the analysis 
procedure that assigns an energy above threshold to D^D^n^ or D^D^'j events below the 
threshold [13]. The quantum numbers l^"*" imply that the X(3872) has an S-wave coupling 
to D*^D^. Its tiny energy relative to the D*^D^ threshold implies that it is a resonant 
coupling. This system is therefore governed by the universal behavior of particles with short- 
range interactions and an S-wave threshold resonance that is predicted by nonrelativistic 
quantum mechanics [2] . The universal properties of the system are determined by the pair 
scattering length a only. If a > 0, one of the universal properties is the binding energy: 
Ex = l/(2/ia^), where fi is the reduced mass of the pair. Another universal property is the 
root-mean-square separation of the constituents: rx = a/\/2. Identifying the mass difference 
in Eq. (2) as —Ex, we find that the charm mesons in the X(3872) have an astonishingly 
large rms separation: rx = 4.9^}^4^ fm. 

The universal properties of an S-wave threshold resonance in the 2-body sector can be 
derived from the universal transition amplitude for scattering of the constituents, which is 
a function of the total energy E of the pair relative to the scattering threshold in their 
center-of-momentum frame: 

A(E) = J^/^ . (3) 

— 1/a -I- V — 2/i£/ — ze 



If a is positive, this amplitude has a pole at the energy of the bound state: E = —Ex, 
where Ex = l/(2yua^). The rate for a process involving the bound state can be calculated 
diagrammatically by introducing a vertex for the coupling of the bound state to its con- 
stituents. Up to an irrelevant phase, the coupling constant G for the vertex is the square 
root of the residue of the pole in the amplitude in Eq. (3) ai E = —Ex'- 

G=[M'Ex/^i'Y'\ (4) 

Taking into account the amplitude 1/v^ for the constituents of X to be D*'^D'^ or D'^D*^ 
from Eq. (1), the vertex factor for the coupling of X to either pair of charm mesons is G /\/2. 



III. CHARM MESON SCATTERING INTO Xtv 

The masses of many different particles enter into the cross section that we will calculate. 
It is therefore useful to introduce a compact notation for these masses. We denote the masses 
of D*^ and D*^ by M^q and M=^+, or by M* if the difference can be neglected. We denote 
the masses of D^ and D^ by Mq and M_|_, or by M if the difference can be neglected. We 
denote the masses of vr'' and tt"'" by mo and m+, or by m,r if the difference can be neglected. 
The mass of the X(3872) is Mx = M*o + Mq — Ex, where Ex is its binding energy. The 
differences between the D* mass and the sum of the D and tt masses are denoted by 

(5) 
(6) 

We will refer to the energy scale set by these mass differences as 5. 

In XEFT, there is an extremely small energy scale: the binding energy of the X(3872), 
which is less than 1 MeV. In the generalization of XEFT that includes charged charm 
mesons and charged pions, there is also the relatively small energy scale 5. The reason the 
D* mesons have relatively long lifetimes is because 5 is small compared to m,r, which makes 
their hadronic decay rates comparable to their electromagnetic decay rates. A pion with 
momentum comparable to itlt^ is generally referred to as a soft pion. We will refer to a pion 
with kinetic energy comparable to or smaller than 5 as an ultrasoft pion. The momentum 
scale for an ultrasoft pion is \^m^, which is about 30 MeV. 

Two types of small mass ratios that appear in XEFT calculations are the ratio of a pion 
mass to a charm meson mass, such as itio/Mq = 0.072 and tjiq/M^q = 0.067, and the ratio 
of a mass difference 5 to a pion mass, such as 5o+/mo = 0.043 and Soo/mQ = 0.053. The 
numerical values of these ratios are all about 5%. We therefore treat all the ratios m,r/^(*) 
and 6/m.,r as order e. The ratio of the binding energy of the X(3872) to its mass is much 
smaller: Ex/Mx = 0.0001 for Ex = 0.42 MeV. We will treat this ratio as 0{e^). We 
will organize our calculations to include all terms suppressed by 0(e), but we will feel free 
to drop corrections of O(e^). Thus a multiplicative factor of the mass of X(3872) can be 
approximated as Mx ~ M^o + Mq. The reduced mass of D*'^ and Z)° can be expressed as 
H ?» MM^,/Mx, since the corrections are 0{Ex/Mx) and hence O(e^) in our power counting. 
Likewise, the reduced mass of 7r+ and X is fi^^+x ~ Mxiti^/{2M^) and the reduced mass of 
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and D^ is fi-j^+o ~ Mm^/M^,. The corrections to these approximations are suppressed 
by 6/M(^^) and hence O(e^). 

We would like to calculate the cross section for the scattering of an ultrasoft pion and the 
X(3872) into a pair of spin-1 charm mesons D*D*. We find it convenient to first calculate 



FIG. 1: Feynnian diagram for the reaction D*~^D*^ — )■ Xir^ . Spin-1 and spin-0 charm mesons are 
represented by sohd+dashed and sohd hnes, respectively. Pions are represented by dashed hnes. 
The X(3872) is represented by a sohd+dashed+sohd hne. 

the cross section for the time- reversed process: the formation of X(3872) by the reaction 
D*D* — )■ Xn. At leading order in the pion transitions, this reaction proceeds only through 
the Feynman diagram in Fig. 1. Our notation for the lines representing the X, the charm 
mesons, and the pion is such that the number of solid lines and dashed lines is conserved 
at every vertex. The pion is ultrasoft if the collision energy of the D* and D*, which is 
their kinetic energy in the center-of-mass frame, is small compared to rrij^. In the diagram 
in Fig. 1, the virtual D is off its energy shell by an amount that is small compared to 771^^. 
Thus this reaction is within the region of validity of XEFT. The reactions D*D — )■ Xtt 
and DD — )■ Xtt can proceed through diagrams similar to the one in Fig. 1 except that the 
number of dashed lines is not conserved at every vertex. The virtual charm meson that 
is exchanged must therefore be off its energy shell by an amount of order 171^^. Thus these 
reactions can not be described by XEFT. 

We proceed to calculate the cross section for the reaction /)*+/)*o — ). Xn^. We take the 
incoming momenta of the D*^ and D*^ to be ±p, and we take the outgoing momenta of the 
7r^ and X to be ±fc. We take the collision energy Ecm = p'^/M^, to be small compared to 
m^, so the vr"^ is ultrasoft. The T-matrix element is 



"TT) 



iG/V2)ig/U)ik.e+)ie*^.eo) 
Ex+[p-{MjMx)k]y{2fi)' ^' 

where e*x, e+, and eo are the polarization vectors of the X, D*^, and Z)*°, respectively. 
The D*Dn coupling constant g/ f-n- can be determined from the measured partial width for 
D*+ -^ D%+: 

g/U = 2.8 X 10"^ MeV-^/l (8) 

The coupling constant g defined here differs from the dimensionless coupling constant g = 0.6 
used in Refs. [3, 4] by a factor of y/2m^. The simple expression in the denominator of 
Eq. (7) for the propagator of the virtual D° can be obtained most easily by exploiting the 
approximate Galilean invariance of this reaction. Galilean invariance holds when the sum 
of the masses is exactly equal in the initial and final states. For the reaction D*+l)*^ — ). 
Xtt^, the sum of the masses decreases between the initial and final states by only about 
(5o+ ~ 5.85 MeV, which is about 1 part in 700. We can take advantage of the approximate 
Galilean invariance by calculating the D propagator in the rest frame of the X(3872). The 
momentum p — {M^/Mx)k in the denominator of Eq. (7) is just the momentum transferred 
through the virtual D meson in that frame. The propagator calculated in the center-of- 
mass frame involves the sum or difference of three kinetic energies. Those three terms can 
be combined to give the single kinetic energy in the denominator of Eq. (7) plus terms that 
have been neglected because they are suppressed by a factor of Ex/Mx = O(e^). 



Squaring the T-matrix element in Eq. (7) and averaging and summing over the spins of 
the D*+, D*^, and X, we get 

1 V- 1^,2 ^ 2G\g/Urfi'k^ 

9^^J 3[2fiEx+p'-2{M,/Mx){p-k) + {M,/Mxyk^f' 



After integrating over the momenta of the outgoing X and tt, our final result for the reaction 

rate is 

2G\g/U)^fi''fi^xk^ 

37rA(p2,(M2/M|.)A;2,-2/iEx) 

where A is the triangle function: A{a,b,c) = a'^ + b'^ + c^ ~ 2{ab + be + ca). The pion 
momentum k is determined by conservation of energy: 






^' 5,^ + Ex + ^. (11) 



By completing the square in the variable p^ in the denominator of Eq. (10) and dividing by 
the relative velocity v^ei = 2p/M=^, we obtain our final expression for the cross section: 

a\D*+D*^ . Xtt+I = G\g/UfMxMlm^k^/{2A'np) 

^ J \p^-{Mj2M){mJ,^-MxEx)f + MlmjQ+MxEx/M^- ^ ' 

At small momentum j9, the cross section diverges as 1/p but the reaction rate v^c\<^[D*^ D*^ — )■ 
Xtt^] is well-behaved. The asymptotic behavior of the cross section for large momentum p 
is 



,*o ^ ^^+1 _^ {g/i-fMxM{mJ^,f l\2MxExf'' 

12p2 



a[D*+D*' ^ X7r+] -^ yyu^l^^^x^^^y..^,^, y....x^x, ^^3^ 



In Fig. 2, we show the reaction rate ■Urcio'[-D*+Z)*° — )> Xtt+J as a function of the D*^D*^ 
collision energy E^.^ = p'^/M^, for several values of the binding energy Ex of the X(3872). 
From the expression in Eq. (12), we can see that there are two competing momentum 
scales that govern the behavior of the reaction rate at small p: (m,r'5o+)^^^ ~ 30 MeV and 
{MxExY^'^, which is approximately 60 MeV if Ex = 1 MeV. If the binding energy Ex 
decreases to less than about ijIt^Sq^/Mx ~ 0.2 MeV, the peak in the reaction rate shifts 
from zero collision energy to a positive value near (m^(5o+ — MxEx)/{2n). This is the 
collision energy for which there is no momentum transferred through the virtual D^. The 
peak at nonzero collision energy is evident in the curve for Ex = 0.125 MeV in Fig. 2. 

The reaction rate for B*^D*^ — )■ Xtt" can be calculated in a similar way. In addition 
to the Feynman diagram in Fig. 1, which involves exchange of a virtual D°, there is a 
second diagram that involves exchange of a virtual D^ . The T-matrix element is therefore 
the sum of two terms. The term corresponding to the Feynman diagram in Fig. 1 differs 
from the expression in Eq. (7) by an isospin Clebsch-Gordan factor of l/\/2- In the energy 
conservation condition in Eq. (11), (5o+ must be replaced by 5qq. 

We now consider the reaction tt^X — )> D*^D*^. We take the incoming momenta of the 
7r+ and X to be A; and we take the outgoing momenta of the D*^ and D*^ to be p. They 
are related by the conservation of energy condition in Eq. (11). The collision energy is 
-^cm = k'^/i'^fJ'wx)- The energy threshold for the breakup process is 5o+ + Ex- The cross 
section can be obtained from the production cross section in Eq. (12) by changing the flux 
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FIG. 2: Reaction rate v^e\a for D*+]J*^ — ). Xi:^ as a function of the collision energy E^ra = p^ /M^, 
for several values of the binding energy Ex ■ 



factor from M*/(2p) to ii-wx/k and by changing the phase space factor from /x^x^/tt to 
M,p/(27r): 



a[XTi^ 



*+ n*oi 



D*^D 



G\g/UfMxMlm^kp/{2A7i) 



[p2 - (M,/2M)(m^5o+ - MxEx)f + MlmJo+MxEx/M^' 



(14) 



IV. -kX elastic scattering 

Another process in the DDtttt sector that is calculable using XEFT is tt^X elastic scat- 
tering. At leading order in the pion transitions, vr+X elastic scattering proceeds through 
six one-loop diagrams. The two diagrams in Fig. 3 involve virtual D*^ and D*^ mesons 
and virtual D^ and D~ mesons, respectively. There are also four additional diagrams with 
Weinberg- Tomozawa vertices. One of them is the first diagram in Fig. 3 with the virtual 
D*^ propagator shrunk to a point to obtain a vr+Z)" — 7r+D° vertex. Another one is ob- 
tained from the second diagram by shrinking the virtual D^ propagator to a point to obtain 
a Ti^D*^ — 7r+Z)*° vertex. The other two diagrams involve a tt^D*^ — n^D*^ vertex and a 
Ti^D^ — Ti^D^ vertex, respectively. 

In Ref. [3], the authors carried out a power-counting analysis for XEFT that determined 
that the D*—Dn vertices can be treated perturbatively. They did not consider the Weinberg- 
Tomozawa (WT) vertices for Dir — Dtt and D*tt — D*7r. These come from the kinetic term 
in the HHxPT Lagrangian, 



C = Ti[Hi{iDo)baHt 



(15) 



FIG. 3: Two of the six one- loop Feynman diagrams for tt^X elastic scattering. The other four 
diagrams involve a Weinberg- Tomazawa vertex and add up to zero. 

where the chiral covariant derivative is 

(L'o)6a = d^dba - ^k> dQ-n^a + 0(7r*^) , (16) 

and our definitions for the fields iff,, H\, and tx can be found in Appendix A of Ref. [3]. 
The WT coupling for D^ mesons is obtained by plugging the second term in Eq. (16) into 
Eq. (15), yielding 

^WT = ^47r-^^7r+Po, (17) 

where Pq is the field for the D*^ . There is a similar coupling for the D*^. Next we rewrite 
Eq. (17) in terms of nonrelativistic pion fields, as is appropriate for the XEFT Lagrangian. 
The relativistic fields tt^ when written in terms of nonrelativistic fields are 



TT 



± _ _^^ (e-^™-*7r± + e*™-*7r^1") , (18) 



where the tt^ (^^^) denotes a nonrelativistic field that annihilates (creates) tt^ mesons. The 
XEFT Lagrangian for the WT vertex is 

C-wT = ^Po^(7r+t7r+ - 7r-t7r-)Po , (19) 

with a similar term for D^* mesons. We have kept only the terms in which the phase factors 
^±tmTyt Qg^^cel. Other terms describe processes that are outside the range of validity of XEFT. 
The WT vertex is 0{Q^) in the power counting of Ref. [3]. This means a diagram with a WT 
vertex is the same order as a diagram obtained by replacing the WT vertex with a virtual D 
or D* propagator and two D* — Dn couplings, since D* — Dtt vertices are 0{Q) and D and 
D* propagators are 0{1/Q'^). The one-loop diagrams for 7r"''X(3872) elastic scattering with 
a WT vertex are then the same order in Q as the diagrams in Fig. 3. However, applying 
charge conjugation to Eq. (19), we see that the Feynman rule for n^D'^^*^ — )■ tx^D^^*'> has 
the opposite sign as the Feynman rule for -k^D^^*' — )■ ti^D^^\ so the four diagrams with 
a WT vertex add up to zero. The only nonvanishing one-loop contribution to 7r^X(3872) 
elastic scattering comes from the diagrams in Fig. 3. 



An explicit one-loop calculation of tt~^X elastic scattering would be worthwhile as a test 
of the calculational technology of XEFT. For collision energies below the break-up threshold 
(5o+ + Ex, the scattering amplitude will be real- valued. As the binding energy Ex decreases 
to 0, the constituents of the X(3872) have a larger and larger probability of being well 
separated. The cross section near threshold should therefore reduce in this limit to the sum 
of the cross sections for scattering off the individual constituents. 

V. SUMMARY 

XEFT is an effective field theory that was originally designed to describe systems con- 
sisting of D*^D^, D^D*^, and D^D'^n^ with total energy relative to the L)*°D° threshold 
that is small compared to the 8 MeV D*^D^ threshold. The only important energy scale 
in this effective field theory is the binding energy Ex of the X(3872). XEFT can be gen- 
eralized to an effective field theory that includes charged charm mesons and charged pions 
and describes systems consisting of D*D, DD*, and DDtt with total energy relative to the 
D*D threshold that is small compared to m.^. In addition to the tiny energy scale Ex, this 
effective field theory also describes the ultrasoft energy scale 6 set by the difference between 
D* — D mass splittings and 771,^. We have pointed out that XEFT can also be applied to 
systems consisting of D*D*, D*D7r, DD*tt, and DDtttt with total energy relative to the 
D*D* threshold that is small compared to 171,^. We have used XEFT to calculate the cross 
section for tt+X — )> /}*+/)*o for ultrasoft collision energy. This cross section is completely 
determined by the masses and widths of the charm mesons, the pion masses, and the binding 
energy of the X (3872). 
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